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Abstract: The logarithmic broadening predicted by the systematic low-energy effective 
field theory for the confining string has recently been verified in numerical simulations of 
(2 + l)-d SU{2) lattice Yang-Mills theory at zero temperature. The same effective theory 
predicts linear broadening of the string at low non-zero temperature. In this paper, we 
verify this prediction by comparison with very precise Monte Carlo data. The comparison 
involves no additional adjustable parameters, because the low-energy constants of the ef- 
fective theory have already been fixed at zero temperature. It yields very good agreement 
between the underlying Yang-Mills theory and the effective string theory. 
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1. Introduction 

At low temperature, in Yang-Mills theory static quarks and anti-quarks are confined to 
one another by an unbreakable string. The string itself is an interesting object with pe- 
culiar dynamical properties. The relevant energy scale for these dynamics is set by the 
string tension a, which defines the asymptotic slope of the linearly rising quark-anti-quark 
potential. The string tension also controls the amplitude of transverse fluctuations of the 
string. As was pointed out a long time ago by Liischer, Symanzik, and Weisz [1], these 
fluctuations are described by a systematic low-energy effective string theory. In fact, the 
transverse fluctuations are massless modes (i.e. Goldstone excitations) of a spontaneously 
broken translational symmetry of the string world-sheet, also known as capillary waves in 
condensed matter physics. The effective string theory plays the role of chiral perturbation 
theory for these Goldstone modes. Using the effective theory, Liischer has calculated the 
leading correction to the linear quark-anti-quark potential 



where d is the space-time dimension and (d — 2) counts the directions transverse to the 
string world-sheet [1, 2]. Accurate numerical simulations have shown the validity of that 
expectation [3-12]. The effective theory also accounts for the string width [13]. Due to 
its transverse fluctuations, the string broadens as the quark sources are separated. The 
cross-sectional area swept out by the fluctuating string increases logarithmically with the 
distance r between the static sources, i.e. 



Here wi {r/2) is the leading-order string width at the distance r/2, half-way between the 
external static quark sources, and ro is a length scale that enters the effective string theory 
as a low-energy parameter. The above expression for the width is universal and applies to 
strings in Yang-Mills theory as well as to fluctuating interfaces between different phases of 
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condensed matter. It has been verified in numerical simulations for the interfaces separating 
the two low-temperature phases of the 3-d Ising model above the roughening transition [14]. 
Similar results have been obtained in a (2 + l)-d Z(4) gauge theory [15]. Since numerical 
lattice calculations of the string width in non-Abelian gauge theories are computationally- 
very challenging, early calculations were affected by large statistical errors and did thus 
not yield conclusive results [16]. Only recently, using the very efficient Liischer-Weisz 
multi- level algorithm [6, 17], the logarithmically divergent string width has been accurately 
verified by simulating (2 + l)-d SU{2) lattice Yang-Mills theory [18]. Recent results on the 
effective string description of the color flux tube in (2 + l)-d SU(N) Yang-Mills theory can 
be found in [19, 20]. 

At low temperature T = 1/(3, the extent (3 of the periodic Euclidean time direction 
enters as another length scale in the low-energy effective string theory. At low T <C -</a, 
the temperature effects on the string tension and the string broadening are again accessible 
to a systematic treatment. At leading order, the temperature-dependent string tension is 
given by 

Tr(d-2) 

Similarly, the leading-order string width at a distance r/2>/5 takes the form [21, 22] 

Ml >/2) = ^.og^ + ^, + (e-^). (1.4) 

Hence, at non-zero temperature the confining string broadens linearly with the distance 
between the external static quarks. Numerical evidence for linear string broadening at finite 
temperature was presented for fluctuating interfaces in the 3-d Ising model in [21, 22] and, 
recently, for (3 + l)-d SU{2>) Yang-Mills theory in [23]. In this paper, using the Liischer- 
Weisz multi-level simulation technique, we perform very accurate Monte Carlo calculations 
of the temperature-dependent string width in (2 + l)-d SU(2) Yang-Mills theory. The 
numerical results are then compared with the corresponding analytic predictions of the 
systematic low-energy effective string theory at the 2-loop level [24]. Since the low-energy 
parameters of the effective theory were previously fixed by a comparison with Monte Carlo 
zero temperature, there are no adjustable parameters. We find perfect agreement 
between the numerical data and the analytic predictions, which lends further support to 
the quantitative validity of the effective theory. 

The rest of this paper is organized as follows. Section 2 describes the analytic calcula- 
tions in the low-energy effective string theory, which are then compared with the numerical 
results obtained in section 3. Section 4 contains our conclusions. 

2. Analytic Results of the Low-Energy Effective String Theory 

Let us consider the systematic low-energy effective theory describing the transverse fluc- 
tuations of the string world sheet in (2 + l)-d Yang-Mills theory. In terms of the height 
variable h(x,t) the corresponding Euclidean action is given by 



S[h] = [ dt [ 
Jo Jo 



dx°- 



d^hd^h - i (d^hd^h) 2 
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Here x G [0, r] and t £ [0, /?] parameterize the 2-d base-space. Since the string ends at static 
quark charges separated by a distance r, its fluctuation field obeys the boundary condition 
/i(0, t) = h(r, t) = 0. The boundary conditions explicitly break translation invariance and 
one might have expected boundary terms to be present in the effective action as well [6]. 
Remarkably, due to open-closed string duality, such terms are absent and the prefactor of 
the first sub-leading term is uniquely determined [25]. This behavior even extends to the 
next order six-derivative term. Indeed, up to that order the effective action coincides with 
the one of the Nambu-Goto string [26-28]. 

The confining string of a lattice Yang-Mills theory suffers from lattice artifacts which 
disappear only in the continuum limit. Away from the continuum limit, in the strong cou- 
pling regime, the lattice string is rigid, i.e. it follows the discrete lattice steps and does not 
even support massless fluctuations. When one crosses the roughening transition at interme- 
diate values of the coupling, one enters the rough phase with massless Goldstone excitations 
of the string, to which the low-energy effective string theory applies. Since the lattice the- 
ory is invariant only under discrete rotations and not under the full Poincare group, before 
one reaches the continuum limit two additional terms proportional to 5^u=i 2 d^d^hd^d^h 
and to Ylw=i 2^^) i en t er the effective theory in the bulk. Since these terms contain four 
derivatives, they are of sub-leading order. Hence, they have no effect on the Liischer term 
or on the leading behavior of the string width. As a result, the Liischer term is indeed com- 
pletely universal. Provided its world-sheet is rough, even a lattice string supports exactly 
massless modes which contribute — 7r/24r to the static quark potential. In particular, this 
term is not affected by lattice spacing artifacts, since such effects are of sub-leading order. 
In the following, we will work sufficiently close to the continuum limit of the lattice theory. 
In that case, the effective action of eq.( |2.1| ) without the additional non-rotational-invariant 
terms is sufficient. 

Recently, an analytic expression for the string width at zero temperature has been 
derived from the effective string theory at the 2-loop level [24] 
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The parameter q = e 



2-KlT 



depends on the ratio r = z/3/(2r) and 





is the first Eisenstein series. The leading-order expression for the width 
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contains the Dedekind 77-function 

00 

V(r)=q^l\(l-q n )- (2-6) 

71=1 

This function and the Eisenstein series E 2 obey the following transformations under the 
inversion r — > — r _1 

V(r) = -X=V(~) , Mr) = ~ 2 E 2 (--) - — . (2.7) 

y —IT T T T VKT 

Using these relations, the expression of eq.( |2.2D in the zero-temperature limit, /3 S> r, is 
readily converted into an expression in the limit r ^> f3. In this case, the temperature is 
still small, but the string length r is much larger than the inverse temperature j3. 

For completeness, we present also the expression for the force F at finite temperature, 
defined as the first derivative of the potential. At leading order the force is given by 

f-(r)-' + ^-^~)- ("J 

The potential at next-to-leading order has been obtained from the computation of the 
partition function using the effective string action of eq. (|2.l| ) [29] . Then the force is given 
by the expression 



F(r) = F lo (r) 



IT 2 d 



r 2E 4 (--) - E 2 2 {--) - + -Ae 2 (--) 
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72o-/3 4 dr 

This equation shows that the 2-loop string tension at finite temperature is given by 



(2- 
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3. Numerical Results for (2 + l)-d SU(2) Lattice Yang-Mills Theory 

Since it is least problematical in numerical simulations, we concentrate on (2 + l)-d SU (2) 
Yang-Mills theory on a cubic lattice of size L x L x /3, with the Euclidean time extent j3 
determining the inverse temperature. The partition function is then given by 

Z= IvU exp(-5[C/]) =T\ ! dU Xtlt exp(-S[U]). (3.1) 

J Z ,JSU(2) 



X,fJ, 



Here dU x> u denotes the local gauge invariant Haar measure for a parallel transporter vari- 
able U x> n £ SU(2) located on the link (x,/x). We use the standard Wilson plaquette action 

S W =4E nU^U x+ ^Ul +p ,,Ul u ], (3.2) 

where g is the bare gauge coupling and fi is a unit-vector pointing in the /i-direction. As 
usual in lattice calculations, all physical quantities are expressed in units of the lattice 
spacing which is set to 1. 
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The external static quarks at the two ends of the confining string are represented by 
Polyakov loops 



IK 



x+t2,2 



(3.3) 



which wind around the periodic Euclidean time direction. We choose the 2-direction to 
represent Euclidean time. When two external quarks are separated in the 1-direction, the 
world-sheet of the confining string extends in the space-time directions [i = 1,2, while the 
string fluctuates in the transverse 3-direction. The static quark potential V(r) is extracted 
from the Polyakov loop correlation function 

($ $ r ) = | Jw $ $r exp(-S[U]) ~ exp(-pV(r)), (3.4) 

in the zero-temperature limit f3 — > oo. 

The force between the two static quarks is extracted from the discrete derivative of 
the static potential 

F(f) = V(r + 1) - V(r) = -- log , (3.5) 

where, writing explicitly the lattice spacing a, f = r + | + 0(a 2 ) is an improved definition 
of the point where the force is computed [30]. 

The width of the string is determined from the connected correlation function 

C{X3) = l$^T ~ <Pa> ' (3,6) 

of a pair of Polyakov loops with a single plaquette 

P x = Tr[U x ,iU x+ij2 Ul +iA Ul 2 ], (3.7) 

in the 1-2 plane, which measures the color-electric field along the 1-direction of the string 
as a function of the transverse displacement x%. The plaquette is located at the site 
x = (r/2,t,xs) and thus measures the transverse fluctuations of the color flux tube at the 
maximal distance r/2 from the external quark charges. 

We have performed numerical simulations at the bare gauge coupling A/g 2 = 9.0 on a 
90 2 x 16 lattice, corresponding to T/T c 0.38. We have used the Liischer-Weisz technique 
[6, 17] with a single level and with slices of thickness 4. In table || we report the number 
of sub-lattice updates for various ranges of distances between the static sources. 

In figure [l], we plot the force between the two static quarks as a function of the 
distance between the sources. The dotted line is the force in the leading order approx- 
imation of eq. ( |2.8D . The solid line is the next-to- leading order approximation given by 
eq. ( |2.9D . It is important to note that there is no free parameter and a = 0.025897(15) is 
the zero-temperature string tension. The next-to-leading order expression already provides 
a very good description of the numerical data. Considering the accuracy of the data, it is 
reasonable to expect that the residual tiny discrepancy can be described by the next-to- 
next-to-leading order correction. 



- 5 - 



r range 


updates 


4-8 


3200 


9-15 


24000 


16-22 


32000 


23-25 


80000 



Table 1: Number of sub-lattice updates of the single level Liischer-Weisz algorithm for various 
ranges of distances between the static sources. 
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Figure 1: The force F(r) as a function of the distance between the two static quarks. The dotted 



line and the solid line are the leading order approximation of eq.(2.8) and the next-to-leading order 



approximation of eq.(2.9), respectively. 



For the measurement of the string width, the correlation function is fit to the ansatz 



A exp(-x§/s) 



1 + B exp(— x\/s) 



+ K, 



(3.8) 



($o$ r > " v d/ 1 + D exp(-x§/s) 

and the squared width of the string is then obtained as the second moment of the correlation 



function 



w 2 {r/2) 



f dx 3 x\ C(x 3 ) 



(3.9) 



/ dx 3 C(x 3 ) 

In figure we show the bell-shaped correlation function ($>Q$> r P x ) / ($o^r) as a func- 
tion of the transverse displacement of the plaquette from the plane of the string world-sheet. 
The data refer to a pair of static sources at fixed distance r = 17. The solid line is a fit of 
the data using the ansatz ( pljj) . 

The scale that determines the logarithmic broadening of the string at zero temperature 
was previously determined as ro = 2.26(2) such that r^y/a = 0.364(3) [18]. It is important 
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Figure 2: The ratio ^$^> P ^ as a function of the transverse displacement £3 of the plaquette from 
the plane of the string world-sheet. The distance between the two static sources is r = 17. The 
solid curve is a fit using eq.(^). 

to note that a and r§ are the only low-energy parameters entering the 2-loop formula for 
the width of the flux tube. Hence the theoretical prediction resulting from the low-energy 
effective theory is completely determined without any additional adjustable parameter. 

In figure ^ we report the observed dependence of the flux tube width on the distance 
between the static sources as resulting from numerical Monte Carlo simulations. The two 
solid lines represent the fully constrained prediction of the low-energy effective string theory 
corresponding to the choices tq = 2.26 ±0.02. The numerical data show the expected linear 
behavior for sufficiently large separations of the sources and the agreement with the analytic 
formula turns out to be excellent. 

Finally, it is important to notice that the next-to-leading expression ( |2.2j ) for the width 
is not sensitive to the next-to-leading correction to the string tension. In fact, the 2-loop 
result for the width goes up to corrections of order l/(cr/3 2 ) and so the corrections of order 
1/(<t/3 2 ) 2 to the string tension will appear as 3-loop corrections to the width. Looking at 
the data for the width in figure ||, the accuracy is such that we are not very sensitive to the 
next-to-leading correction to the string tension and one cannot detect the tiny difference 
that the direct measurement of the force in figure [j] shows. 

4. Conclusions 

Using the Liischer-Weisz multi-level simulation technique, we have performed high-accuracy 
numerical simulations of the confining string connecting two static quarks in (2+l)-d SU (2) 
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Figure 3: The squared width of the confining string at its midpoint w 2 {r/2) as a function of the 
separation of the static sources r. The dotted line is the leading order prediction of the low-energy 
effective string theory using ro = 2.26. The two solid lines represent the 2-loop prediction using 
r = 2.26 ±0.02. 

Yang-Mills theory. The results for the finite temperature string width have been compared 
with analytic 2-loop calculations performed in the corresponding systematic low-energy 
effective string theory. Since the low-energy parameters a and tq were already determined 
at zero temperature, there were no further adjustable parameters in the comparison of the 
effective theory with the numerical simulations of the lattice Yang-Mills theory. Remark- 
ably, there is perfect agreement between the underlying Yang-Mills theory and the effective 
string theory, thus confirming the quantitative correctness of the latter. With sufficient 
numerical resources (which are still modest on today's standards in lattice field theory) it 
is certainly feasible to extend our study to other gauge groups and to (3 + 1) dimensions. A 
similar investigation in (3 + l)-d SU(3) Yang-Mills theory would definitely be possible. In 
particular, it would be interesting to investigate whether the effective theory works equally 
well in the numerically accessible regime as in the (2 + l)-d case investigated here. 
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